Kolmogorov and Kelvin- Wave Cascades of Superfiuid Turbulence at T = 0: 

What is in Between ? 
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As long as vorticity quantization remains irrelevant for the long-wave physics, superfiuid turbu- 
lence supports a regime macroscopically identical to the Kolmogorov cascade of a normal liquid. At 
high enough wavenumbers, the energy flux in the wavelength space is carried by individual Kelvin- 
wave cascades on separate vortex lines. We analyze the transformation of the Kolmogorov cascade 
into the Kelvin- wave cascade, revealing a chain of three distinct intermediate cascades, supported by 
local-induction motion of the vortex lines, and distinguished by specific reconnection mechanisms. 
The most prominent qualitative feature predicted is unavoidable production of vortex rings of the 
size of the order of inter-vortex distance. 

PACS numbers: 67.40.Vs, 47.32. Cc, 47.37.+q, 03.75.Kk 
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Nowadays, superfiuid turbulence [l], 0] — a structured 
or non-structured tangle of quantized vortex lines — is 
attracting much attention stimulated, in particular, 
by advances in experimental techniques allowing studies 
of different turbulent regimes in diverse superfiuid sys- 
tems, such as 4 He @, Sj, 3 He-B [a, 0, and Bose-Einstein 
condensates of ultacold atoms @, 0|- ^ n superfiuids at 
T = 0, vorticity can only exist in the form of topolog- 
ical defects — vortex lines of microscopic thickness, the 
circulation of velocity around which being equal to the 
liquid-specific quantum k. Speaking generally, the dy- 
namical mechanisms governing superfiuid turbulence are 
fundamentally different from those of classical turbulence 
(see, e.g., recent review Q and references therein). 

A new wave of interest in dynamics of superfiuid 
turbulence came with the experiment by Maurer and 
Tabeling |8|, who observed that superfiuid turbulence 
in 4 He formed by counter-rotating discs is indistinguish- 
able from classical turbulence at large length scales, in 
particular, exhibiting the classical Kolmogorov cascade. 
Shortly, the same effect was found in superfiuid turbu- 
lence generated by a towed grid l 9|. In the experiments 
[1, @] , the fraction of normal component is considerable 
making analysis of vortex tangle dynamics and structure 
significantly complicated [lfj [Hj]. (Considerations re- 
garding possible energy spectra in this case are presented 
in Ref. [12] •) However, the similarity between classical 
and superfiuid turbulence exists even at practically zero 
temperature, which was first observed in numerical sim- 
ulations [3, [3 EB| , and, just recently.for the first time 
confirmed by measurements in 3 He-B 6]. 

By the nature of a cascade regime, implying that the 
kinetic times get progressively shorter down the hierar- 
chy of length scales, instantaneous structure of turbu- 
lence follows the evolution at the largest length scales 
(typically of order of system size), where the energy flux 
(per unit mass) e is formed. At very low temperatures, 
due to the absence of frictional dissipation, the flux e 
must be carried down to scales significantly smaller than 
the (related to e) typical separation between the vortex 
lines Iq. At small enough length scales, the energy flux 
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FIG. 1: Spectrum of Kelvin waves in the quantized regime. 
The inertial range consists of a chain of cascades driven by 
different mechanisms: (1) reconnections of vortex- line bun- 
dles, (2) reconnections between nearest-neighbor vortex lines 
in a bundle, (3) self-reconnections on single vortex lines, (4) 
non-linear dynamics of single vortex lines without reconnec- 
tions. The regimes (3) and (4) are familiar in the context of 
non-structured vortex tangle decay 0, EH • 



is carried by pure Kelvin-wave cascades on separate vor- 
tex line s lla . [TtJ , the cutoff being due to sound radiation 
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The fact that superfiuid turbulence at large compared 
to Iq length scales may be consistent with the classical 
Kolmogorov law is not surprising (a formal proof is men- 
tioned below). It is well known [l| that macroscopic ve- 
locity profile of a rapidly rotated superfiuid mimics solid- 
body rotation, which is accomplished by formation of 
a dense array of vortex lines aligned along the rotation 
axis. By the same mechanism, "stirring" a superfiuid 
one can produce vorticity in the course-grained up to 
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length scales larger than lg superfluid velocity field, in- 
distinguishable from that of a normal fluid, the underly- 
ing vortex tangle being organized in polarized "bundles" 
of vortex lines. What turns out to be a puzzle [20], how- 
ever, is how the vortex tangle looks like when one zooms 
in down to scales of order of the interline separation Iq, 
where the vorticity is essentially discrete. 

In this Letter, we analyze the structure of turbulence 
at all length scales tracing the transformation of the clas- 
sical regime, described by the Kolmogorov law at large 
length scales, into the quantized regime, in which the dis- 
creteness of vortex lines is important, in the fundamental 
case of zero temperature. The analysis relies on the large 
parameter 

A = ln(Z /ao) > 1, (1) 

where ao is the vortex core radius. In realistic He ex- 
periments, A ~ 15. The attention to the problem of link- 
ing the two regimes was drawn recently by L'vov et al. 
[20], who realized that it is impossible to directly cross 
over from Kolmogorov regime to the pure Kelvin-wave 
cascade, and put forward the idea of a bottleneck, with 
specific dynamical implications. The Achilles' heel of the 
treatment of Ref. [2(| is taking it for granted [l(| E| that 
the coarse-grained macroscopic description of quantized 
vorticity remains valid down to the scale of Iq. 

We show that the locally-induced motion of the vortex 
lines radically changes dynamical picture already at the 
scale 

ro ~ A^lo , (2) 

with the interline separation related to the energy flux 
by 

lo ~ (A/c 3 /e) 1/4 . (3) 

In the range of wavelength ro > A > A* , 

A* = 2oM 1/2 , (4) 

there takes place a chain of three cascade regimes, in 
which the energy flux e is carried by locally-induced mo- 
tion combined with vortex line reconnections. The three 
regimes are distinguished by their specific types of re- 
connections: (i) reconnections of vortex-line bundles, (ii) 
reconnections between nearest-neighbor vortex lines, (iii) 
self-reconnections on single vortex lines — the mechanism 
introduced earlier by one of us [2l| in the context of 
the decay of non-structured superfluid turbulence. The 
existence of the regime (iii) means unavoidable produc- 
tion of vortex rings of the typical size A* at a rate im- 
mediately following from ([3]) by conservation of energy. 
Namely, ~ kA 1 / 2 //^ rings are emitted per unit time in 
the characteristic volume of Iq. Remarkably, this rate is 
~ A 3 / 2 times smaller than the rate of vortex ring pro- 
duction characteristic of non-structured superfluid tur- 
bulence psm. 



For realistic values of A, sharp distinction between the 
three sub-regimes is likely to be lost, although character- 
istic features of strong turbulence, such as generation of 
a spectrum of vortex rings by the mechanism (iii), might 
manifest themselves. At the wavelength scale A*, self- 
reconnections cease and the weak-turbulent regime sets 
in, with purely non- linear Kelvin- wave cascade [lq . This 
regime covers a significa nt p art of the inertial range until 
eventually, at the scale 

A ph ^[A 27 («/c) 25 ^] 1/31 (5) 

(c is the sound velocity) , the cascade is cut off due to the 
radiation of sound by Kelvin waves. With k/c~ do we 
have A p h «< A,. 

The structure of turbulence is summarized in Fig. |T]). 
We emphasize that the notion of energy spectrum E(k), 
where E(k)dk gives the energy per unit mass associated 
with variations of fluid velocity over length scales ~ fc _1 
in the interval dk, is practically meaningful only in the 
classical regime. In the quantized regime, the relevant 
degrees of freedom are waves on vortex lines, while even 
a perfectly straight vortex line has a nontrivial spectrum 
E{k) ~ k. On the experimental side, recently introduced 
vortex line visualization technique in 4 He [j] could pro- 
vide the most direct probe for the quantized regime. 

Analysis. The key to understanding vortex dynamics 
at zero temperature is given by Kelvin- Helmholtz's theo- 
rem, which states that a vortex line element moves with 
local fluid velocity. Mathematically, this is reflected in 
the Biot-Savart equation [l[, 

s = v(s), v(r) = £- J (s - r) x ds /|s - r| 3 . (6) 

Here v(r) is the superfluid velocity field, s is the time- 
evolving radius- vector of the vortex line element, the dot 
denotes differentiation with respect to time, the vector 
so has the same physical meaning as s, understood as 
an integration variable, and the integration is along all 
the vortex lines. In the long-wave limit, when discrete- 
ness of vortex lines is irrelevant, Eqs. (JB]) can be coarse 
grained: The second equation turns into Biot-Savart re- 
solvent restoring the field v from its curl, w. The first 
equation coarse-grains into Euler's equation for incom- 
pressible fluid, 

W a = -VffdpWa + WpdpVa . (7) 

Here the first term in the r.h.s. comes from translation of 
the vortex array by the flow, the second term originates 
from bending of the array by velocity gradient. This 
formally proves the equivalence of (structured) superfluid 
and normal-ideal-incompressible-fluid turbulence. For a 
non-structured tangle, coarse graining trivially leads to 
v = 0. 

For our purposes, it is instructive to formally decom- 
pose the integral © into the self- induced part, v SI (s), 
for which the integration is restricted to the vortex line 
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containing the element s, and the remaining contribution 
induced by all the other lines v x (s), 



v(s) =v SI (s)+v I (s). 



(8) 



Since the velocities v SI and v 1 define the r.h.s. of Eq. ©, 
the competition between them is crucial for the problem. 

The leading contribution to v SI is given by the local in- 
duction approximation (LIA) , which reduces the integral 
over the vortex line to its local differential characteristics, 



v SI ( s ) = A^f s' x s", 



A R = HR/a ) , (9) 



where the prime denotes differentiation with respect to 
the arc length and R is the typical curvature radius. A 
necessary condition of applicability of the LIA is A^, 3> 1. 
Taking into account that A^ is a very weak function of 
R, we shall treat it as a constant of typical value Ar ~ A. 
Note, however, that despite the fact that the condition 
(JTJ) is typically well satisfied, using the LIA is not always 
appropriate. — Being an integrable model, LIA does not 
capture reconnection-free (purely non-linear) kinetics of 
Kelvin waves [l6| • 

To determine the crossover scale ro, consider the struc- 
ture of the vortex tangle in the classical regime. By the 
definition of ro, at length scales r > ^ turbulence mim- 
ics classical vorticity by taking on the form of a dense 
coherently moving array of vortex lines bent at curva- 
ture radius of order r. Velocity field of this configuration 
obeys the Kolmogorov law 
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(10) 



Here and below the subscript r means typical variation 
of a field over distance ~ r. On the other hand, the 
value of v r is fixed by the quantization of velocity circu- 
lation around a contour of radius r, namely v r r ~ Kn r r 2 , 
where n r is the areal density of vortex lines responsible 
for vorticity at the scale r. Note, that scale invariance 
requires that on top of vorticity at the scale r there be a 
fine structure of vortex bundles of smaller sizes, so that, 
mathematically, n r r 2 is the difference between huge num- 
bers of vortex lines crossing the area of the contour r in 
opposite directions. The quantity n r is related to the flux 
by 
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r > r . 



(11) 



The underlying dynamics of a single vortex line in the 
bundle is governed by v l r and vf 1 . While, by its definition, 
v\, ~ v r , which is given by Eq. (|10|) . the self- induced part 
is determined by the curvature radius r of the vortex line 
according to Eq. ©, 
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At length scales where v\ 3> vf 1 , the vortex lines in the 
bundle move coherently with the same velocity ~ uL 



However, at the scale ro ~ (A 3 k 3 /e) 1 / 4 , the self-induced 
motion of the vortex line becomes comparable to the col- 
lective motion, vf l ~ v\. At this scale, individual vortex 
lines start to behave independently from each other and 
thus ro gives the lower cutoff of the inertial region of the 
Kolmogorov spectrum iflQ]) . 

Since ro is the size of the smallest classical eddies, the 
areal density of vortex lines at this scale is given by the 
typical interline separation, n ro ~ 1/Zq. With Eq. I|ll[). 
we arrive at jU)-®. 

At the scale ro, turbulence consists of randomly ori- 
ented vortex line bundles of size ro, left by the classical 
regime. The typical number of vortices in the bundle 
is given by n ro rQ ~ A. The length ro plays the role of a 
correlation radius in the sense that relative orientation of 
two vortex lines becomes uncorrelated only if they are a 
distance > ro apart. On the other hand, the crossover to 
the quantized regime means that each line starts moving 
according to its geometric shape, as described by Eq. ((9]). 
Therefore, reconnections, at least between separate bun- 
dles, are inevitable and, as we show below, capable of 
sustaining the flux e. 

Reconnections play the leading role at ro > A > A*. 
Although this region is relatively narrow as compared to 
the whole Kelvin- wave inertial range, it is significantly 
large in the absolute units. Before going into the details 
of the reconnection-assisted regimes, we describe the re- 
maining and dominant region of the cascade. As was 
shown by the authors [l6j |. at a sufficiently small wave- 
length, a strongly turbulent cascade of Kelvin waves is 
replaced by a purely non-linear cascade, in which the re- 
connections are exponentially suppressed. The spectrum 
of Kelvin- wave amplitudes bk, k ~ A , in the non-linear 
cascade has the form 



b k = (e/K 3 p) 1 /io fe -6/s ! 



(13) 



where O is the flux of energy per unit vortex line length 
supported by the non-linear cascade. The value of A* 
can be determined by matching the energy flux e with 



Q/pll, where b k 



A*. With Eq. ©, we then 



obtain Eq. g}. 

Kelvin waves decay with emitting phonons [lOj . For 
Kelvin waves of wavenumber ~ k, the pow er of sound 
emission per unit line length is given by 19] 



n« 



A 6 K s pbik n /c 5 



(14) 



This dissipation mechanism is negligibly weak all the way 
down to wavelengths of order A p h, given by Eq. |5]), 
where Hk/plg becomes comparable to e. The scale 
Aph A* gives the lower dissipative cutoff of the Kelvin 
wave cascade. 

Now we focus on the strongly turbulent regimes at 
ro > fc _1 > A*. The key quantity here is the energy 
transferred to a lower scale after one reconnection of vor- 
tex lines at the scale A; -1 , which, following Ref. [2l|], can 
be written as 



£fe - /(7)A/)K 2 fc- 



(15) 
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Here, /(■y) is a dimensionless function of the angle 7 at 
which the vortex lines cross, 7 = corresponding to par- 
allel lines. Its asymptotic form is 

/( 7 )~ 7 2 , 7 «1. (16) 

Although, at the scale r , there is already no coupling 
between vortex lines to stabilize the bundles, they should 
still move coherently — the geometry of neighboring lines 
at this scale is essentially the same over distances ~ ro — 
until the whole bundles cross each other. It is possible, 
however, that vortex lines within the bundle reconnect. 
One can show that such processes can not lead to any 
significant redistribution of energy, and thus to a defor- 
mation of the bundle at the scale ro, because they happen 
at small angles so that the energy (fT5|) is too small. In- 
deed, the dimensional upper bound on the rate at which 
two lines at distance I <C ro can cross each other is, from 
Eq. ([9]), Ak/tqI, while the actual value should be much 
smaller due to the strong correlations between line ge- 
ometries. Taking into account that the number of lines 
in the bundle is (to/Iq) 2 and that 7 ~ //ro, the contribu- 
tion to the energy flux from these processes is bounded 
by (l/ro)e. 

Crossing of the bundles results in reconnections be- 
tween their vortex lines and Kelvin waves with somewhat 
smaller wavelength A are generated. The coherence of the 
initial bundles implies that the waves on different vortex 
lines must be generated coherently. Thus, at the scale 
A < ro, vortex lines should be also organized in bun- 
dles of length A that are bent with the amplitude of the 
Kelvin waves bk, k ~ A , while the correlation radius 
for vortex line configurations in the transversal direction 
is ~ bk. Then, reconnections between the bundles at the 
scale A transport the energy to a lower scale. The cascade 
of bundles should repeat itself self-similarly in a range of 
wavelength Iq <C k ,bk < ro, in which the notion of 
bundles is meaningful. The spectrum of Kelvin waves bk 
in this regime can be obtained from the condition ik = £, 
where ik is the energy flux per unit mass transported by 
the reconnections at the scale given by 

ik ~ (k/pbl) N k e k t-\ (17) 

Here, we take into account that the correlation volume is 



b 2 /k, the Nk ~ (bk/lo) 2 is the number of vortex lines in 
the bundle, and t^ 1 ~ Auk 2 is the rate at which the bun- 
dles cross. Physically, bk determines the typical crossing 
angle, 7 ~ bkk, thereby controlling the energy lost in one 
reconnection. Thus, the spectrum of Kelvin waves in the 
bundles is 

b k ~ r^k- 2 . (18) 

At the wavelength ~ Ab = A 1 ' 4 Iq, the amplitudes be- 
come of order of the interline separation bk ~ Iq and the 
cascade of bundles is cut off. At this scale, bkk <C 1, 
so that the mechanism of self-reconnections is strongly 
suppressed. The kinetic times of the purely non-linear 
regime are too long to carry the flux e [16| . We thus con- 
clude that at A c < A < At the cascade is supported by 
nearest- neighbor reconnections, the amplitudes bk being 
defined by the condition of constant energy flux per unit 
length and the crossover scale A c being associated with 
the condition bk^\/\ r ~ A c meaning that at A < A c the 
self-crossing regime takes over. The observation crucial 
for understanding the particular mechanism of the cas- 
cade and thus finding bk is that each nearest-neighbor 
reconnection (happening at the rate cx A/ A 2 per each 
line element of the length ~ Ah) performs a sort of parallel 
processing of the cascade for each of the wavelength scales 
of the range [A c , Ab]. For the given wavelength scale 
A ~ 1/k, the energy transferred by a single collision is 
oc A(bkk) 2 X, and with the above estimate of the collision 
rate per the length Ah, this readily yields the estimate 
bk ~ io(A&fe) -1 ' 2 , and, correspondingly A c ~ Iq/A 1 / 4 . 

In the range A c 3> k^ 1 ^ A*, the cascade is driven 
by self-reconnections of vortex lines giving the spectrum 
bk ~ k~ l [2l| . This regime is replaced by the purely 
non-linear regime in the vicinity of ~ A* (the actual 
transition region may be rather wide [16l|). 

To conclude, the transformation of classical-fluid Kol- 
mogorov cascade of superfluid turbulence into the pure 
Kelvin-wave cascade requires three intermediate stages 
associated with locally-induced motion and reconnec- 
tions of vortex lines, as illustrated in Fig.[TJ 
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